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1.0 INTRODUCTION 


The NASA Lewis Research Center employs a general computer program (Reference 1) 
for calculating the thermodynamic performance of jet propulsion engines. To cal- 
culate off-design engine performance, the user must input component maps. These maps 
define the characteristics of the various components over their full range of 
operating conditions . 

For advanced propulsion systems these characteristics are not generally known. 
Furthermore, the typical user of the program is not sufficiently knowledgeable and/or 
cannot afford the time to do an extensive design analysis of the component in question. 
Instead he usually scales some available map. 

The objective of the study is an improved method of representing the turbine 
component when performing calculations of off-design performance for advanced air- 
breathing jet engines. This method, which is a computer program called PART, is com- 
patible in both form and format with the cycle program of Reference 1 and the example 
map representation of Reference 2. 

The current program is a follow-on to NASA Contract NAS3-21999. Under the 
original contract an axial flow turbine model for large flow size machines was de- 
veloped. Under the current contract, the model was extended to include both small 
axial flow turbines (i.e., flow sizes down to about 1 pps) , and small fixed geometry 
uncooled radial flow turbines. 

Because this report contains a description of the input-output data, values of 
typical inputs, and sample cases, it is suitable as a user’s manual. A brief 
description of the engineering analysis used to generate the program is given near 
the end of the report. 

The program uses turbine design point data as input to generate off-design 
values of turbine flow-function and total-to-total efficiency over a range of 
pressure ratios and speeds specified by the user. A user-specified option will also 
permit calculating design point cooling flows for the axial flow machines, and the 
corresponding change in turbine efficiency. The cooling flow subroutine, developed at 
the Lewis Research Center, is described in Reference 3. 
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Table I. Summary of Turbine Design Point Data. 


No. 

Data Base Name* 

No. of Stages 

DHQTD ** 

tffd** 

i 

HPT 1-1 

i 

0.0596 

14.6 

2 

HPT 1-2 

i 

0.0705 

16.4 

3 

HPT1-3 

i 

0.0335 

88.5 

4 

HPT2-4 

2 

0.0670 

17.3 

5 

HPT2-5 

2 

0.0787 

32.4 

6 

HPT 3-6 

3 

0.0810 

45.5 

7 

LPT 1-1 

1 

0.0220 

45.0 

8 

LPT 1-2 

1 

0.0425 

45.0 

9 

LPT2-3 

2 

0.0571 

58.5 

10 

LPT 2 -4 

2 

0.065 

60.4 

11 

LPT4-5 

4 

0.0665 

106.0 

12 

LPT4-6 

4 

0.0709 

134.4 

13 

LPT 6- 7 

6 

0.0814 

104.9 

14 

PT3-1 

3 

0.0800 

210.0 

15 

AT3-1 

3 

0.0590 

— 

16 

AT3-2 

3 

0.0785 

— 

17 

AT3-3 

3 

0.0635 

43.16 

18 

AT4-4 

4 

0.0499 

38.85 

19 

VAT1-1 

1 

0.044 

99.0 

20 

VAT1-2 

1 

0.060 

60.0 

21 

VAT1-3 

1 

0.0238 

290.0 

22 

VAT1-4P 

1 

0.0328 

61.8 

23 

VAT1-4X 

1 

0.0328 

61.8 

24 

VAT2-5P 

2 

0.0636 

61.8 

25 

VAT2-5X 

2 

0.0636 

61.8 

*HPT 

- High Pressure Turbine 



LPT 

- Low Pressure Turbine 



PT 

- Power Turbine 




AT 

- Air Turbine Test Rig 



VAT 

- Variable Area Turbine 




** Symbols defined in Table III. 











Table II. Summary of Variable Geometry Turbines 
Included on Data Base. 


Turb ine 
No. 

Designat ion 

No. of 
Stages 

First Stage Nozzle Area Ratios 
Z 

i 

VATl-l 

1 

50.0 62.5 75.0 87.5 100.0 

2 

VAT1-2 

1 

71.0 86.0 100.0 109.0 120.0 

3 

VAT1-3 

1 

76.0 84.0 92.0 100.0 108.0 116.0 

4 

VAT1-4P 

1 

70.0 100.0 130.0 

5 

VAT1-4X 

1 

70.0 100.0 130.0 

6 

VAT2-5P 

2 

70.0 100.0 130.0 

7 

VAT2-5X 

2 

100.0 


3 













The Aircraft Engine Group of the General Electric Company has a turbine 
data base consisting of 25 turbines having design point turbine flow functions 
ranging from about 14 to 290. The number of stages for each of these turbines 
together with the approximate design point values of specific work output divided 
by inlet total temperature (DHQTD) and flow function (TFFD) are summarized in 
Table 1. The last sevel turbines shown in the table are variable-geometry turbines. 
Table II shows the set of first stage nozzle area ratios for each of the seven 
variable-geometry turbines. Five of these variable-geometry turbines were generated 
by turbine design and off-design computer programs similar, if not identical, to 
that described in Reference 4. Two of the turbines shown in Table II were generated 
from air turbine test carried out by the Lewis Research Center. The results of 
these tests are given in References 5 to 8. The Table II designation of the NASA 
test turbines have been given a trailing X. The analytical prediction of the perfor- 
mance of these turbines obtained from Reference 9 has been a trailing P. 

All of the turbines in data base discussed above are of the large axial-flow 
type. In order to obtain data for small axial-flow turbines and radial flow turbines 
the open literature was used. References 13 thru 17 give the references used for 
small axial flow turbines. References 18 thru 24 give those for radial turbines. 
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2 . 0 PROGRAM STRUCTURE 


A flow chart showing the flow of control in the NASA parametric turbine 
program is shown in Figure 1, After the input has been read and processed , 
the program carries out a simple pitch line analysis starting with the last 
stage of the turbine. The analysis starts at the exit of the turbine stage 
(in order to avoid iteration) and calculates the bucket and nozzle flow 
angles. This stage geometry is then used to generate the stage flow and loss 
characteristics using the analytically based correlations developed during 
the program. Successive stages are then calculated until the first stage is 
reached. The first stage characteristics are then generated, and the stages 
stacked for each value of the first stage turbine nozzle area specified. If 
the turbine is cooled, then the procedure given in Reference 3 is used to 
calculate both the cooling flow requirements, and the coolmd turbine effi- 
ciency. Finally, the output is processed to obtain a turbine map representa- 
tion compatible with the cycle deck of Reference 1. 
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Figure 1. Flov Chart Showing Flow of Control in Parametric 
Turbine Program. 












3 . 0 PROGRAM INPUTS 


or.’io'?.'/! * 

OF POO;? 


rj 


i * 


y 


All of the PART input* ere of the free-field format (NAMELIST) type, end 
begin in column two. There it no specified order to the inputs. The program 
initially lists the contents of the NAMELIST INPUT together with the default 
settings of all the input variables. The user may then change as many of the 
inputs at desired. The program then echoes the updated NAMELIST. If none of 
the inputs are changed, the program will execute the first example case and 
the user can inspect the output. The input variables together with the de- 
fault settings for a turbine are summarized in Table III. 

The first six input variables in Table III are used to control the number 
and values of speed, pressure ratio, and nozzle area ratios (transformed into 
nozzle angle) to be written on the output files. For example, in the input to 
the first example case shown in Figure 2, all of the corrected speed and pres* 
sure ratio arrays are used by the program, but only the first three positions 
in the area ratio array. Note that speeds and pressure ratios are entered in 
increasing value, but that area ratios are entered in decreasing value (this 
is so that the nozzle angles will be written in increasing order on the output 
file). Speeds less than 102 should not be used. The input to the second ex* 
ample case shown in Figure 3 illustrates the use of the first six variables 
to limit the size of the output files. 

Some of the design point inputs will be calculated internally by the pro* 
grso, if the user inputs the correct value to trigger the calculation. This 
subset of inputs together with the required settings are swmmarized in Table 
IV. An input value equal to or less than zero will trigger all of the calcu- 
lations with the exception of exit swirl angle, here a value greater than 90 
degrees must be input (180 degrees is recommended). 

A minimum set of design point input would consist of NSTG and DHQTD. 

The values of TFFD and XNRTD could be input as 100.0 and the resulting values 
of TFFD interpreted as percent. The user could then use the settings in Table 
IV to trigger program calculations of the remaining design point information. 
The use of the program to calculate the number of stages will frequently re* 
suit in a single-stage turbine, since the only upper limit on turbine rsdius 
is the limiting value of rim speed. This is not usually sufficiently restric* 
tive to require the use of additional stages. 

If the user wishes the program to calculate the value of design point 
cooling flow, and the corresponding decrease in turbine efficiency, the JCOCL 
switch in the NAMELIST INPUT should be set to 1. The program will then list 
the contents of the NAMELIST INPUT 1 together with the default settings of all 
the variables. The user may then change as many of the inputs as desired. 
Since the default settings of the NAMELIST INPUT are for an uncooled turbine, 
these inputs (namely, TON, PTIN) must be changed in order to successfully 
calculate cooling flows. The input variables for NAMELIST INPUT1 together 
with their default settings are siznmarized in Table V. The input to the 
second example case shown in Figure 3 illustrates the proper format for a 
cooled three-stage turbine. 
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Table til. Default Betting* for Variable* in Namelist "Input." 



Units 

Default 

Values 

(Axial Turbine 

^ Description 

RSPDS 

Rone 

15 

Number of Speed Lines Desired 

APCNC(l) 

Rone 

I0Z to 150X 

The Array of Percent Corrected Speeds (Max of 15) 

NPR 

Rone 

20 

Rusher of Pressure Ratios Desired 

APR(l) 

Rone 

1.1 to 4.6 

The Array of Pressure Ratios (Ha* of 20) 

RAR 

Rone 

3 

Number of Firat Stage Nozzle Area Ratios 

ARN 

Rone 

1.3, 1.0, 0.7 

Array of Rottle Area Ratios (Ha* of 6) 

HSTG 

Rone 

1 

Number of Turbine Stages (Ha* of 6) 

JCOOL 

Rone 

0 

Cooling Flow Switch (0*Uncooled; l*Cooied) 

DHQTD 

Btu/lba *R 

0.03270 

Specific Work Output Divided by Inlet Temperature 

ETATTD 

Rone 

0.923 

Turbine Total-to-Total Efficiency 

TFFO 

Ibm *R 1/2 
nec*psi 

62.98 

Turbine Inlet Flow Funct ion f i.e., 
(TFF - W/i7/P t ) 

XRRTD 

rpm/'R 1 ^ 2 

193.52 

Turbine Corrected Speed (XNRT • 

PSIO 

Rone 

0.8511 

Average Turbine pitch line loading, l*e. f 

FSID - <DHQT/(2(U//r7) 2 /g 0 Jl/NSTG 

« 

ANCSWX 

Degrees 

15.2 

Em it Pitch Line Swirl Angle (Positive When Opposite 
to Direction of Rotation) 

XMZXD 

Rone 

0.373 

Exit Fitch Line Axial Mach Number 

TTIH 

• R 

518.67 

Turbine Inlet Total Temperature 

rriH 

ps is 

14.696 

Turbine Inlet Total Pressure 

FARGD 

Rone 

0.0 

Turbine Inlet Fuel-Air Ratio 

ETAN 

None 

.94 

Nozzle Efficiency (Ratio of Exit Actual to ideal 
Kinetic Energy) 

R3QR2 

None 

1.0 

Pitch Line Radius Ratio (Rotor Exit to Inlet) 

RHQRT3 

None 

.733 

Radius Ratio at Rotor Exit (First Guess) 
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■PARAMETRIC TURBIHE' 


P “A*R*T 


AXIAL OR RADIAL TURBINE? 

-M 


1 -RADIAL, 0 OR CR-AXIAL 


DESCRIPTION OF INPUT VARIABLES IN PART PROCRAM 
( NAML LIST INPUT) 

ITS : TTIN-DEC R , PT 1 N-PSIA ,W-PPS , H-BTU/LBH , W-RPH 

MAP range t ar ra y values must be in ascending order) 


NSPDS 

APCNC 

NPR 

APR 

NAR 

ARN 

NSTG 


* i 5 


■20 


NO 0 f SPEEDLINES DESIRED! HAX* I 5 ) 

AKRAf Of PERCENT CORRECTED SPEEDS 
NO Of PRESSURE RATIOS 0 £ S I RE D ( MAX- 20 ) 
ARRAV Of PRESSURE RATIOS 
NO Of NOZZLE AREA RAT I OS ( MAX*6 ) 

ARRAY Of FIRST STAGE NOZZLE AREA RATIOS 
NO Of TURBINE STAGE S < MAR *6 > 


INTEGER SWITCHES 

JCOOL -0 CO ClING fLOW SWITCHI0-UNCOOLED, I-COOLED ) 


DESIGN 
OHQTD 
ETATTD 
Jf f D 

KNRTD 

PSID 

ANGSWX 

xmZ xO 

TT1N 

PTIN 

FARGO 

ETan 


POINT VA h, UE S Of s 

SPECIFIC WORK OUTPUT DIVIDED BY TT1N 
TURblNl TOTAL-TG-TuTAL EFFICIENCY 
TURBINE INLET KlOW FUNCTION 
4 TF F -W-SuRl ( TT1N >/PTlN ) 

TUKblNt COKRECVED SPEEDl N/SGRTI TTIN ) ) 
AVERAGE TURBINE P 1 TCH LINE LOAOINC 
PS1G*DmGT/( 2 * ( U/SURTiTTlN ) • *2 /G J ) I/NSTG 
EXIT PITCH LINE SWIRL ANGLE ( ♦COUNT-ROT ) 
EXIT PITCH LINE AXIAL MACH NUMBER 
TURBINE INLET TOTAL TEMPERATURE 
TURttlNL INLET TOTAL PRESSRUE 
T UR b 1 Nt INLET FUEL-AIR RATIO 
NOZZLE EFFiRATIO ACTUAL TO 10EAL EXIT KE ) 


GEOMETRY 
R30RZ 
RHOR 1 3 


SPECIF ICAtlONS : 

P1TCHLINL RAO RAT J OI ROTOR EXIT To INLET) 
RAOiUS RATIO AT ROTOR E X I T< F 1RST GUESS) 


NAME L 1ST 


INPUT 


NSPDS 

APCNC 

1 

5 

9 

13 

NPR 

APR 

1 

5 

9 


• IS, 

( I >* 

1* .0JUJ . 

S4 .Jdkikixi , 

9 4 . add* , 
I 34 .hi*** , 

* Z4, 

< I )* 

1 . 1400 . 
) . 7004 . 
2 . 4 000 . 


24.4404, 
66 . 0404 , 

l 00.0000, 
1 40 . 0000 , 


i . 2000 . 

i . bduu , 

Z . 68*0 , 


30.0000, 
70.4404, 
1 10.4000, 
1 S0 . 0400 , 


I . 4000 , 
2 .0000 » 
2.8000, 


40 

B0 

120 


Figure 2. Input to First Example Case 


0000 , 
0000 , 
. 0000 , 


. 6000 , 
. 2000 , 
. 0000 , 



ORIGINAL PASS r C: 
0F p OOR QUALITY 


13 3.2000, 3.4000, 

l? 40000. 4.20U0* 

NAR - 3, 

ARN ( 1 1 0 MM 

1 1.3000. 1 . 0000 . 

5 0. 0. . 

1 1.2490. 1.1440. 

NSTG * 1 , JCOOL ■ 0 . 

OHQTO • 0.032780, 

jpfD a 62.980000* 

p <, | Q m 0.851100, 

ANGSWX* IS • 2 00000 . 

p j 1 h a 1 4 . to 9 6000 , 

tlAN a 45.940000, 

RHGRT3- 0 733J00. 

£ HU NAME LIST INPUT 

ENTER CnANuES TO NAMELIST INPUT 
•SlhPUTS 

NAMELIST INPUT 

NSPDS » 15, 

APCNC t 1 )■= 

1 10 


ETATTO- 
XNRTD » 
KMZXD - 
TTIN - 
FARGD • 
R3QR2 - 


0000 , 
0000 , 
0000 
0000 


5 50 . 

9 90 . 

13 130 

HPR - 20, 

APR < I >* 

1 1 . 1000 , 

b 1 . 7000 , 

9 2.4000. 

13 3.2000, 

17 4.0000, 

NAR - 3 . 

ARN III* 

1 1.3000, 

5 0 . 

1 1.2490, 

NSTG - 1 . JCOGL • 

DHUTD * 0. 032700, 

q . o r . 9b 0000 , 

P510 • 0.8S1100. 

ANGSUX- IS. 200080. 

P J 1 a 14.690000, 

£ TAN • 0.940000, 

IRnQRT3« 0.733300, 

£ NO NAME LIST INPUT 

NASA OUTPUT ON IfC-15. 16 


20 . 0000 , 
60 . 0000 , 
100.0000 , 
140.0000, 


. 2000 , 
. 8000 , 
, 6000 . 
. 4000 , 


4.2000. 


I .0000 , 

0 . 

1 . 1 440, 

0, 


ETATTD 

XNRTD 

XMZXD 

TTIN 

F ARGD 

R3QR2 


3,6000. 3. 

4.4000. *• 

0.7000. *• 

1.0000, *‘ 


0.923000 
193 . 520000 

0. 373000 
518. 669990 

0 . 

1 . f 00000 


30.0000. 

40 

70.0000. 

80 

110.0000, 

1S0.0000. 

120 

1.4000. 
2.0000. 
2.8000. 
3 . 6000 . 

4 . 4000 . 

1 

2 

3 

3 

4 

0.7000, 

4 

1 .0000, 

1 


0.923000. 

193.520000. 

0.373000, 

510.669998, 

0 . 

1 . 000000 , 


8000, 

6000, 


0000 . 


. 0000 . 

. 0000 . 

. 0000 . 


. 6000 . 
.2 000. 
. 0000 , 
.8000. 
, 6000 . 


. 0000 , 


Figure 2 


(Continued) Input to First Example Case 



OPr'M r ’. 

Cr 


PARAMETRIC TURBINE 


P*A*R*T 

AXIAL OR RADIAL TURBINE? 1* RADIAL , » OR CR-AXIAL 

•0 

DESCRIPTION OF INPUT VARIABLES IN PART PROGRAM 
(NAMELIST INPUT > 

UNITS:TTIN-DEG R .PT1N-PSIA.W-PPS .H-BTU/LBM.N-RPM 

MAP RANGE ( ARRAY VALUES MUST BE IN ASCENDING OROCR ) 
NSPDS -IS 
APCNC 
NPR -20 
APR 

NAR *6 * 

ARM 

NSTC -6 


NO OF SPEEDLINES DE S I RE D< MAX- 1 5 ) 

ARRAY OF PERCENT CORRECTED SPEEDS 
NO OF PRESSURE RATIOS DE S I R£D< MAX-20 > 
ARRAY OF PRESSURE RATIOS 
NO OF NOZZLE AREA RAT I OS < MAX-6 > 

ARRAY OF FIRST STAGE NOZZLE AREA RATIOS 
NO OF TURBINE STAGE SI MAX-6 ) 


INTEGER SWITCHES 

OCOOL *0 COOLING FLOW SWITCH! 


0-UNCOOL E D, 1 -COOL ED 


) 


DESIGN 

DHQTD 

ETATTD 

TFFD 

XNRTD 

PSIO 

ANGSWX 

XMZXO 

TTIN 

PTIN 

FARGO 

ETAN 


POINT VALUES OF i 

SPECIFIC WORK OUTPUT DIVIOED BY TTIN 
TURBINE TOTAL -TO-TOTAL EFFICIENCY 
TURBINE INLET FLOW FUNCTION 
< TFF-W*SQRT( TTIN >/PTIN ) 

TURBINE CORRECTED SPEED! N/ SORT I TTIN ) ) 
AVERAGE TURBINE PITCH LINE LOADING 
PSlO-DHOT/!2*( U/SORT (TTIN)**2/CJ ) >/NSTC 
EXIT PITCH LINE SWIRL ANGLER ♦COUNT -ROT ) 
EXIT PITCH LINE AXIAL MACH NUMBER 
TURBINE INLET TOTAL TEMPERATURE 
TURBINE INLET TOTAL PRESSRUE 
TURBINE INLET FUEL-AIR RATIO 
NOZZLE EFF ( RATIO ACTUAL TO IDEAL EXIT KE ) 


GEOMETRY 

R3QR2 

RHQRT3 


SPECIFICATIONS* 

PITCHLINE RAD RATIO! ROTOR EXIT TO INLET) 
RADIUS RATIO AT ROTOR £XIT(FIRST GUESS) 


NAMELIST INPUT 

NSPDS - 15, 

APCNC < I >- 

1 10 . 0000 , 

5 50.0000, 

9 90.0000, 

13 130. 0000 , 

NPR • 20, 

APR II)- 

1 1 . 1000 , 

5 1 . 7000 , 

9 2.4000, 


20 . 0000 , 
60.0000, 
1 00 . 0000 , 
140.0000, 


1 . 2000 , 
1 . 8000 , 
2.6000, 


30.0000, 

70.0000. 
1 10 . 0000 , 
150.0000, 


1 .4000, 
2 . 0000 , 
2 . S000 , 


40 

80 

120 


1 

2 

3 


0000 , 

0000 , 

0000 , 


.6000. 

. 2000 , 

. 0000 . 


Figure 3. 


Input to Second Example Case 
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3 . 2 000 , 

4.IMI. 

■ 3 . 

I 1 I* 


. 3000 . 


,2490. 


13 
17 
NAR 
ARN 

1 1 

5 * 

1 » 

NSTG - 1 . 

OHCiTD ■ 

TFFO * 

PS1D * 

ANCSUX- 
PTIN - 
ETAN « 

RHQRT3- 
END NAMELIST 
ENTER CHANGES 
•SINPUT 
-NSPDS-6, . 

■APCNC-20. .*40. ,60. .8* 
■APR*2 .0,2. 5, 3. 0,3. 5, 3. 8, 


3,4000. 
4 . 2000 . 


JCOOL - 
0.032780. 
62.980000. 
0.851100. 
15.200000. 
14.696000, 
0.940000, 
0.733300. 
INPUT 

TO NAMELIST 


1 . 0000 . 

0 . 

1 . 1440 , 

0, 


ETATTO" 
XNRTD ■ 
XMZXO • 
TT IN ■ 
FARGO ' 
R3QR2 


INPUT 


,100. ,120- 


• NAR* 1 , 

•ARN* 1 -0, 
-NSTG-3. 
•JCOOL* 1, 
■OHQTD*0 .0635 , 
■ TFF 0*58 . 63 , 
•XNRTD-40.10. 
•PSIO-l .5. 
•ETATTO-0.886. 
•ANGSWX-2.9. 
•XMZXD»0 . 41 , 
•TTIN-2500. . 
•PT 1 N*59 . 8 . 
-FARGU-0.02, 


3.8000, 

4.4000. 


3.8000. 

4.6000. 


0.7000. ■* 

1 . 0000 , 1 . 0000 . 

0. 923000. 
193.520000. 

0.373000, 

518.669998. 

0 . 

1 .000000, 


DESCRIPTION of input variables in part program 

(NAMELIST INPUT!) 

...pof-jN ojoE.t, «« 

tc-total temperature of the cool I NC FLOW 

FARCX-FUEL-A1R RATIO OF THE COOLING FLO MATERIAL 

OF^JS «°EoI Soto. BAADE -AT£«.»A 
(fs^EsKo “?E OF TO.B1HE Al.FO.L 


namelist 
kindof*. 

TC 

TEAR - 
EL1FE - 
END NAMELIST 
ENTER CHANGES 
•SJNPUT1 


INPUT1 
1 10800. 
500-000000. 
I960. . 

0. 1000000E 05 
input i 

TO NAMEL 1ST 


-KINDOF *864000 , 
-TC-700. . 


FARCX 

VEARB 


INPUT1 


namelist input 

NSPOS • 6. 

APCNC t i >- 

1 20 . 0000 . 


40.0000, 


0 . 

1980. . 


60.0000. 


80.0000 • 


Figure 3 - (Continued) Input to Second Example Case 
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5 

1 00.0000, 

120.0000, 

70 . 0000 , 

80.0000 

9 

90.0000, 

100.0000, 

110.0000, 

120.0000 

13 

130.0000, 

140.0000, 

150.0000, 


NPR 

- 5, 




APR 

( I )- 




1 

2.0000, 

2.5000, 

3 . 0000 , 

3 . 5000 

5 

3.8000, 

1 .8000, 

2.0000, 

2 . 2000 

9 

2.4000, 

2.6000, 

2 . 8000 , 

3 . 0000 

13 

3.2000, 

3.4000, 

3 . 6000 , 

3.B000 

17 

4.0000, 

4 . 2000 , 

4.4000, 

4 . 6000 

NAR 

• 1 * 




ARN 

( I )• 




1 

1 .0000, 

1 .0000, 

0 . 7000 , 

0. 

5 

0. 

0 . 



1 

1 .2490, 

1 .1440, 

l .0000. 

1 . 0000 


nstg • 3 , JCOOl - 1, 

DHOTD - 0.063500, ETATTD* 0.086000, 

Tf F D » 50.530000, XNRTD - 40.100000, 

PS1D - • 1.500000, XMZXD • 0.410000, 

ANGSWX- 2.900000, TTIN • 2500.000000, 

PT1N - 59.800000, FARGO • 0.020000, 

ETAN * 0.940000, R3QR2 - 1.000000, 

RHQRT3- 0.733300, 

END NAMELIST INPUT 

PCBLED- 0.06478 PCMCH* 0.02001 EFF4- 0.8835 PRN- 3.297 
NASA OUTPUT ON IFC-15,16 A 


Figure 3. - (Continued) Input to Second Example Case 
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Table IV. Variable Setting* to Trigger Default Calculation 
of Some Design Point Input, 


Variable 

Name 

Setting 

Action Taken 

1 

NSTG 

0.0 

Program Calculates Number of Stages 
(Not Reconnnended) 

ETATTD 

0.0 

Program Calculates Design Point Efficiency 

PS ID 

0.0 

Design Point Loading Set to 0.9 

ANGSWX 

180.0 

Exit Swirl Angle Calculated From Zero Hub 
Reaction 

XMZXD 

0.0 

Sets Exit Axial Mach Number to 0,5 


Table V. Default Settings for Variables in Namelist "Inputl" 


Var iable 
Name 

Units 

Default 

Value 

Description 

KINDOF 

None 

86400000 

An Ordered Combination of Digits 
Representing the Cooling Configura- 
tion of the Turbine 

TC 

* R 

700.0 

Total Temperature of the Cooling 
Flow 

FARCX 

None 

0.0 

Fuel-Air Ratio of the Cooling Flow 

Year 

i 

1980 

First Year of Service for Stator 
Vane Material 

Ye arB 

— 

1980 

First Year of Service for Rotor 
Blade Material 

ELIFE 

hrs 

10,000 

Desired Life of Turbine Airfoil 




c : 


The integer variable RINDOF represents the cooling configuration of the 
turbine. Each blade row starting with the firat stage atator is assigned an 
integer value characterizing the type of cooling employed as follows; 

0 Uncooled 

1 Convection cooling 

2 Convection with coating 

3 Advanced convection 

4 Film with convection (752 trailing edge injection) 

5 Film with convection (502 trailing edge injection) 

6 Film with convection (252 trailing edge injection) 

7 Transpiration with convection (252 trailing edge injection) 

8 Full coverage film 

9 Transpiration 

For example, the 86400000 configuration has the first 'three blade rows 
cooled and the remaining five rows uncooled (a four-stage turbine). For a 
detailed description of the cooling flow calculation and the various cooling 
flow configurations, the reader should consult Reference 3. 
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4 . 0 PROGRAM OUTPUTS 


The basic output from the program consists of two tables. These tables 
show the turbine efficiency and turbine flow function variations for each of the 
first stage nozzle area ratios, pressure ratios, and percent corrected speeds 
specified in the input. The input values of area ratio are converted to first 
stage nozzle angles before being printed out. The output tables for the first 
example case are shown on pages 22 through 27. The table structure is compatible 
with NASA cycle deck requirements given in Reference 2 (pages 23 and 24) . 

The output tables can be visualized as three dimensional, composed of a 
series of planes with each plane assigned a value of nozzle angle, BETA. Then in 
each BETA plane, the dependent variable (ordinate axis) is a function of pressure 
ratio, PR, and corrected speed, rpm. The dependent variables are respectively 
turbine corrected flow, W, and total— to— total efficiency, ETA. 

For example, in the output table on page 30 the forty-five lines of the 
dependent variable correspond to the fifteen values of corrected speed, where each 
speed occupies three lines. And the twenty values of the dependent variable in 
each three line group correspond to the twenty values of pressure ratio. 

In addition to these two tables, there is a terminal listing summarizing 
the results of the cooling flow calculation, if this option was used. The value 
of the total cooling flow, PCBLED, is printed out together with the cooling flow 
for the first stage nozzle alone, PCNCH. The new cooled turbine efficiency value, 
EFF4, is given together with the new value of the total-to-total pressure ratio 
across the turbine, PRN. An example of this printout is shown on page 13 in the 
second example case. With the flows, shaft work, and turbine pressure ratio known, 
the user can calculate the new cooled turbine efficiency, ETATTD, using the book- 
keeping procedure compatible with the cycle deck representation to be employed. 

A cycle deck efficiency scalar could then be used or, if desired, the program could 
be rerun on the uncooled branch using the new design point efficiency value as an 
input . 
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5,0 PROGRAM DIAGNOSTICS 


The PART computer program contains error printouts to aid the user in 
trouble shooting his input. A listing of the error messages and their meanings 
are given below. 

1. LIMITING VALUE OF UTIP=1800.0, CALCULATED VALUE OF UTIP= 

This warning message is printed out if the calculated tip speed of a 
radial turbine exceeds 1800 fps. 

2. LIMITING VALUE OF UHUB=1600.0, CALCULATED VALUE OF UHUB= 

This warning message is printed out only if the calculated rim speed 
exceeds the recommended value (this is a disk stress warning) . 

3. LIMITING VALUE OF ANS=42.0E9, CALCULATED VALUE OF ANS= 

This warning message is printed out if the product of the exit annulus 
area and the rpm squared exceeds the recommended value (this is a centrifugal 
stress limit on the rotor blading). 

4. QUIRE CTR ERROR — (CALLING LINE=,I5,) 

There are eight iterations in the program. Seven of the iterations are 
balanced using the Method of False Position. This method is contained in the sub- 
routine QIREXX. A maximum of 25 passes is allowed for any single iteration to 
balance. If the iteration does not balance within the specified tolerance, the 
error message will appear with the number of the offending iteration in the 15 
Format field. 

Normally, the occurrence of such an error will not cause a problem. 

However, in the case of QIRE loop number five which calculates the turbine 
efficiency for the specified input values of pressure ratio and corrected speed, 
an additional message indicating the convergence error is printed out. This message 
has the form: 

DHQT= , ERR= ,PQP= 

where the blanks contain the current values of specific enthalpy change divided 
by inlet total temperature, the convergence error in pressure ratio, and the 

17 

pressure ratio at which the error occurred. 



The user should inspect the error to see if the degree of convergence is 
satisfactory, if not, it may be necessary to restrict the range of input speeds 
and/or pressure ratios requested. The individual QIRE loops together with the 
calling routines and type of iteration are as follows: 


QIRE 

LOOP 

CALLING 

ROUTINES 

COMMENTS 

1 

INLETX 

Calculates individual stage efficiencies from the input 
value of overall turbine efficiency (NSTG>1) . 

2 

VELRAT 

Obtains the axial velocity ratio across the rotor. 

3 

CHOKEX 

Solves for the value of nozzle Mach number when the 
rotor chokes. 

4 

ROTCKX 

Solves for exit annulus choke location given the 
location of rotor choke. 

5 

PRTEFF 

Calculates efficiency for input values of speed and 
pressure ratio. 

6 

FSTACK 

Solves for the "polytropic" exponent for a multistage 
turbine. 

7 

CHOKEX 

Solves for the value of the nozzle overexpanded Mach 
number after nozzle choke and before rotor choke. 


There is one iteration in the program balanced by the Newton-Raphson Method. A maxi- 
mum of 25 passes is allowed for convergence . If an error occurs the program will 
print out the warning message 

ZERO DERIVATIVE IN NEW RAP LNCALL = I 10 

The value 1 will appear in the 110 field since this is the first Newton-Raphson loop. 


NEWTON- 

RAPHSON CALLING 

LOOP ROUTINE COMMENTS 

1 CHOKEX Calculates the value of UQAT1 at the speed where both 

the nozzle and the rotor are choked. 
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6.0 EXAMPLE CASES 


Two example cases are given in order to illustrate the use of the program. 
The first case utilizes the default settings to generate the output for a single- 
stage, uncooled, variable-geometry turbine. The second case is a single stage radial 
turbine . 

A complete record of the two terminal sessions including a listing of the 
output tables is given on the following pages. The program inputs and outputs have 
been discussed previously in Sections 3.0 and 4.0. 
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ORIGINAL PAG'S ?? 
OF POOR QUALiTs 


•REMOVE CLEARFILES 
•OLD /NASAPT/NASAPART 

•FRN 

*0580 10/20/83 11.122 


‘PARAMETRIC TURBINE’ 


P«A*R*T 

AXIAL OR RADIAL TURBINE? 1-RADIAL, 0 OR CR-AXIAL 

-0 

DESCRIPTION OF INPUT VARIABLES IN PART PROCRAM 
(NAMELIST INPUT) 

UNITS: TT 1N-DEG R . PT I N-P S I A . WPP S . H-BTU/LBM . N-RPM 

«*»«■«»*« 


NSPDS 

APCNC 

NPR 

APR 

NAR 

ARN 

NSTG 


•IS 


•20 


NO OF 'SPEEDLINES DESIRED! MAX-15 ) 

ARRAY OF PERCENT CORRECTED SPEEDS 
NO OF PRESSURE RATIOS DE S I RED< HAX«-0 > 
ARRAY OF PRESSURE RATIOS 
NO OF NOZZLE AREA RAT IOS< MAX-6 ) 

ARRAY OF FIRST STAGE NOZZLE AREA RATIOS 
NO OF TURBINE STAGE S ( MAX-6 ) 


INTEGER SWITCHES 
JCOOL -0 


COOL i NG FLOW SW1TCH(0-UNCOOLED.1-COOLED) 


DESIGN 

DHQTD 

ETATTD 

TFFD 

XNRTD 

PSID 

ANGSWX 

XMZXD 

TT l N 

PTIN 

FARGO 

ETAN 


P ° 1NT SptciFIc°LoRK OUTPUT DIVIDED BY TTIN 
turbine total-to-total efficiency 

TURBINE INLET FLOW FUNCTION 
(TFF*U*$QRT(TTIN>/PT1N> tttn i \ 

TURBINE CORRECTED SPEED< N/SORTITTII N)> 
AVERAGE TURBINE PITCH LINE LOADING 

PSID-OHQT/<2VU/SORT(TTIN)--2/GJ)>/NSTG 

EXIT PITCH LINE SWIRL ANGLE(*COUNT ROT) 

EXIT PITCH LINE AXIAL MACH NUMBER 

TURBINE INLET TOTAL TEMPERATURE 

TURBINE INLET TOTAL PRESSURE 

TURBINE INLET FUEL-AIR RATIO 

NOZZlE EFFIRATIO ACTUAL TO IDEAL EXIT KE ) 


GEOMETRY 

R3QR2 

RHQRT3 


SPI 5!H C Htm"MD RATIO! ^OTOR EXIT TO INLET ) 
RADIUS RATIO AT ROTOR EXlTtFIRST GUESS) 


name list 


INPUT 


NSPOS 

APCNC 

1 

5 

9 

13 

NPR 

APR 

1 

5 

9 


* 15, 

< I >* 

1 0.0*00 , 
50.4004. 
90 . 0000 ' 
130.0000 , 
- 20 , 

( I >« 

1 . 1000 , 
1 . 7000 , 
2.4000, 


20 . 0000 , 
60.0000, 
100 . 0000 , 
140.0000, 


1 . 2000 , 
1 .8000, 
2.6000, 


30 . 0000 . 

70.0000, 
1 10 . 0000 , 
150.0000, 


1 . 4000, 
2 . 0000 , 
2 . 8000 , 


40.0000, 

80.0000, 

120 . 0000 , 


1 .6000, 
2 . 2000 , 
3.0000, 





u r r > 


13 3.20^, 3.4000, 

17 4.0000, 4.2000, 

NAR - 3, 

ARN <l>- 

1 1.3000, 1.0000, 

5 0. 0. , 

1 1.2490, 1.1440, 

NSTG - 1, JCOOL - 0, 

DHQTD - 0.032700, ETATTD- 

TFFD - 62.980000, XNRTO - 

PSID - 0.851100, XMZXD - 

ANGSWX- 15.200000, TT I N - 

PTIN - 14.696000, FARGO - 

ETAN - 0.940000, R3QR2 - 

RHQRT3- 0.733300, 

END NAMELIST INPUT 

ENTER CHANGES TO NAMELIST INPUT 
•SINPUTS 

NAMELIST INPUT 


■ 15, 

( I >* 

10 . 0000 , 

50.0000, 

90.0000. 
130.0000, 

- 20 , 

( I >• 

1 . 1000 . 
1 .7000, 
2.4000. 
3 . 2000 . 
4.0000, 
- 3, 

( I )* 


NSPDS 
APCNC 
1 
5 
9 
13 
NPR 
APR 
l 
5 
9 
13 
17 
NAR 
ARN 

1 1.3000, 

5 0. 

1 1.2490, 

NSTG - 1. 

DHOTD - 
TFFD » 

PSID - 
ANGSWX- 
PTIN - 
ETAN - 
RHQRT 3 * 

END NAMELIST 
NASA OUTPUT ON 


20 . 0000 , 

60.0000, 

100 . 0000 , 

140.0000. 


2000 , 

8000, 

6000, 

4000, 

2000 , 


OCOOL • 
0.032780, 
62.980000, 
0.851 100, 
1 5 . 200000 , 
14.696000, 
0.940000, 
0.733300, 
INPUT 
IFC-15,16 


l .0000, 
0 . 

1 . 1440, 
0 , 


ETATTD- 
XNRTD - 
XMZXD » 
TT 1 N - 
FARGO - 
R3QR2 • 


3.6000, 

4.4000, 

0.7000, 
1 . 0000 , 


3 . 8000 . 
4 . 6000 , 

0 . 

1 . 0000 , 


0.923000, 
193.520000, 
0.373 000, 
516.669998, 

0 . 

1 . 000000 , 


30.0000, 

70.0000, 
110.0000, 
150.0000. 

4 0.0000, 

80.0000, 
120.0000, 

1.4000, 
2.0000, 
2.8000, 
3.6000, 

4.4000, 

1 . 6000 , 
2.2000, 
3.0000, 
3.B000. 
4.6000, 

0.7000, 

0. 

1 .0000, 

1 .0000, 


0.923000, 

193.520060, 

0.373000, 

518.669998, 

0 . 

1 . 000000 , 
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26. 
BETA 3 
RPM 15 
RPM 15 


EFF. 20 
EFF 20 
EFF 20 


EFF 20 
EFF 20 

EFF 

20 

EFF 

20 

F FT " 

“90 


OF POOR Q'JAL 


TURBINE EFFICIENCY VS. 
61 . 68. 

10 . 20 . 

80 90. 




0.2)5 
O. 106 
O. 550 




PR, RPM, AND BETA 
75. 

30. 

100 . 1 




0.209 

0.186 

0.497 




0.204 

0.164 

0.479 


0. 256 
0 . 200 
0.182 
0.464 


0.912 
0. 691 
0 560 


i 

0.750 
0 626 
0. 660 



O. 253 
0 824 
0 723 


0. 879 
0.667 
0 570 


U. JIM 

0 725 
0.617 
0 920 


0 652 
0.905 
0. 785 


0.875 
O. 799 
0 695 


0.807 
O 706 
0.781 


0.649 

0.804 

0.714 


0.643 
0 925 
O 762 


0.925 
0.778 
0. 666 


0.786 
0.697 
0. 695 


0.834 

0.766 

0.706 


0.765 

0.632 

0.553 




0.685 

0.600 

0.878 


0.635 

0.902 

0.743 


0.912 
0. 759 
0.678 


0.768 

0.689 

0.903 


0 . 666 
0.771 
0.696 


0. 746 
0.619 
0 545 




O 667 
O 592 
0 806 


0 626 
0 892 
O. 726 


0. 903 
0.743 
0.670 


0.753 
0.681 
0. 899 


0. 672 
O. 757 
0.691 


0.733 

0.604 

0.538 


0.652 
0.585 
O 852 


0.622 

0.877 

0.712 


0.666 

0.728 

0.663 


0.739 
0.674 
0. 888 


0. 666 
0. 745 
0.664 


0.196 
0 441 


0.371 

0.367 

0.363 

0.360 

0.357 

0 354 


0 677 

0.601 

0. 708 

0.655 

0.637 

0.622 

0 . 598 

0.560 

0.567 

0. 557, 

0.546 

0.534 

0.524 

0.515 


0.711 
0. 591 


0.639 


0 816 


0 84 1 
0. 699 


0. 853 
0.716 


0. 727 


0 858 


0. 845 
0. 734 




0. 730 















jam 


EFF 

20 

0. 690 

0.078 

0.474 

0.422 

0.405 

0. 390 

0 369 

EFF 

20 

0.352 

0.340 

0. 330 

0.321 

0.314 

0. 306 

0 . 303 

EFF 

20 

0 299 

0.295 

0.291 

0.266 

0 285 

0. 282 


EFF 

20 

0.830 

0 731 

0. 62G 

0.569 

0.549 

0.533 

0 . 507 

Err - 

26 

' 0748fi 

VZ ^BB 

h*vt^HB! 

0.451 

0.443 

0 . 4 36 


EFF 

20 

0.424 

0 419 

0.415 

0,411 

0.408 

0.405 


EFF 

20 

0. 893 

0. 627 

0. 736 

0.662 

0. 602 

0. 646 

0 . 620 

EFF 

20 

0. 600 

0 565 

0.573 

0.563 

0 554 

0. 546 

0 . 540 

Bi JiJB 




— <n&? 

0.324 


0. 520 





0.864 

0.614 

0.766 

0. 749 

0.733 

0. 709 




0.676 

0. 663 

0.652 

0 646 

0.641 

0 . 637 




0 630 

O 626 

0.623 

0.619 

IBKmBAmBBI 


miizm 



0.911 



0.813 

|HiWjTajT*V|| 

0.778 

EFF 

20 

0.760 

0. 745 

0.735 

0. 727 

0.721 

0.715 

0 . 709 

EFF 

20 

0.702 

0.696 

0.690 

0.685 

0.680 

0 . 675 


EFF 

20 

0 867 

0.919 

0 900 t 

0.672 

0 859 

0. 849 

0 . 828 

EFF” 


O.S'l T 

0~79!3 

0 709 

0 760 

0.772 



EFF 

20 

0.746 

0.739 

0.733 

0.727 

0. 722 

0.717 


EFF 

20 

0. 817 

0.910 

0.918 

0.901 

0 892 

0.864 

0 . 864 

EFF 

20 

0. 848 

0 636 

0 625 

0.616 

0. 604 

0.794 

0 . 765 



0 "776 





0 . 744 




0.703 

0. 690 

0.925 

0.918 

0.913 

0.906 

0 . 669 



0 874 

0.861 

0.850 

0.838 

0.826 

0.815 

0 805 



0 796 

0.788 

0. 781 

0.775 

0.769 

0.763 




0 676 


■Am 

flH*J£LX±lABii 

0. 924 



EFF 

20 

0.690 

0.877 

0. 665 

0.652 

0 639 

0. 628 

0.818 

EFF 

20 

0.809 

0 801 

0. 794 

0.787 

0.781 

0 . 775 


EFF 

20 

0.089 

0 820 

0.912 

0.927 

0.928 

0.925 

0.913 

EFF" 


oneres 

0 886 

01873 


0.847 



EFF 

20 

0.817 

0.809 

0.801 

0.794 

0.788 

0.782 


EFF 

20 

0.466 

0.774 

0.895 

0.921 

0. 926 

0.925 

0.916 

EFF 

20 

0. 902 

0 889 

0.676 

0. 6G4 

0. 851 

0 840 

0 . 830 


WiBBl 

CTS2A 

BKftlUfli 


0 7758 


0.785 




0.379 

0.721 

0.873 

0.910 

0.919 

0.921 

0.915 



0. 902 

0 689 

0.678 

0.664 

0.852 

0.840 

0.831 



0 822 

0 81 3 

0. 606 

0 799 

0.792 

0. 786 


3 J5S 


0 '.'260 

— cretn — 

0 7845 

O. 695 

0 . 907 

0 . 9 X2 ” 


EFF 

20 

0.897 

0.886 

0.875 

0.661 

0.649 

0.838 

0. 829 

EFF 

20 

0.620 

0 812 

0. 804 

0.797 

0. 791 

0. 785 


EFF 

20 

0 141 

0.596 

0.813 

0.875 

0.891 

0. 900 

0 . 900 

■ ■ me a m 



0“891 


0 870 

07837 




i EFF 

20 

0 616 

0.606 

0.801 

0.794 

0. 786 

0. 782 


5 RPM 

10 

10. 

20. 

30. 

40. 

50. 

60. 

70 . 

£ RPM 

10 

60. 

90. 

100. 

1 10. 

120. 

130. 

1 40 . 

|”RPM” 

io 








1 PR 

20 

1.10 

1 . 20 

* 1.40 

1 . 60 

1 .70 

1 . 60 

2 . 00 

| PR 

20 

2.20 

2.40 

2 60 

2.60 

3.00 

3. 20 

3 . 40 

1 P* 

20 

3.60 

3 80 

4 . 00 

4 20 

4 . 40 

4. 60 
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0.491 
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0.351 

0. 335 

0 323 

0 304 

EFF 


1 T2OT — 

0 "279 

0.270 




0 . 24 f 

EFF 

20 

0.243 
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PARAMETRIC TURBINE 


P-A'R-T 


AXIAL OR RADIAL TURBINE? 1-RADIAL, 0 OR CR-AXIAL 

■ 1 

DESCRIPTION OF INPUT VARIABLES IN PART PROGRAM 
(NAMELIST INPUT) 

UN I TS :TTIN*0EG R , PT I N-P S I A , W-P P S , H-BTU/LBM . N-RPM 


MAP RANGE < ARRAY VALUES MUST BE IN ASCENDING ORDER) 
SPEEDLINES DE S I RE D( MAX- 1 5 > 


PRESSURE RATIOS DE S I RED( MAX-20 ) 
OF PRESSURE RATIOS 
N02ZLE AREA RAT I OS < MAX-6 ) 


JCOOL -0 COOLING FLOW SWI TCH l 0-UNCOOLED , 1 -COOLE D ) 


NSPDS 

• 1 s 

NO OF 

APCNC 


ARRAY 

H PR 

•20 

NO OF 

APR 


ARRAY 

NAft 

■6 

NO OF 

ARN 


ARRAY 

NSTG 

-6 

NO OF 

INTEGER SWITCHES 


DESIGN POINT VALUES OF: 

DHOTO SPECIFIC WORK OUTPUT DIVIDED BV TT1N 

ETATTD TURBINE TOTAL -TO-TOTAL EFFICIENCY 

TFFD TURBINE INLET FLOW FUNCTION 

< TF F -W-SQRT ( TTIN )/PTIN ) 

XNRTD TURBINE CORRECTED S P E E D l N/SQRT( TT I N ) ) 

PS1D AVERAGE TURBINE PITCH LINE LOADING 

PSID-DHOT/(2*l U/ SORT (TTIN )»*2-/GJ ) ) /NSTG 
ANGSUX EXIT PITCH LINE SWIRL ANGLE! *COUNT-ROT ) 

XMZXU EXIT PITCH LINE AXIAL MACH NUMBER 

TTIN TURBINE INLET TOTAL TEMPERATURE 

PTIN TURBINE INLET TOTAL PRESSRUE 

FARGO TURBINE INLET FUEL-AIR RATIO 

£ TAN NOZZLE EFF< RATIO ACTUAL TO IDEAL EXIT KE ) 

GEOMETRY SPECIFICATIONS: 

R3QR2 PITCHLINE RAD RATIO! ROTOR EXIT TO INLET) 

RhORT 3 RADIUS RATIO AT ROTOR EXITtFIRST GUESS) 


NAMELIST INPUT 




NSPDS 

APCNC 

l 

- 15, 

< I >* 

10.0000 , 

20.0000 , 

30.0000, 

40.0000, 

’ 5 

50.0000 , 

50.0000. 

70.0000, 

60.0000, 

9 

90.0000 , 

1 00. 0000. 

1 10.0000, 

120.0000, 

13 

1 30 . 0000 , 

140.0000 , 

150.0000, 


NPR 

APR 

1 

- 20 . 

< I )- 

1 . 1000 * 

1 .2000, 

1 .4000, 

1 .6000, 

5 

1 .7000. 

1 .8000. 

2.0000, 

2.2000, 

9 

2.4 000 . 

2.6000, 

2.8000, 

3 . 0000 , 



13 3.2 m, 

17 4.0040, 

; NAR - 1. 

ARN <I)- 
1 1 • 8888 « 

5 0 . 

1 0.3300 

NSTG - I. 

DHQTD - 
TFFD » 

PSID - 
ANGSWX* 

PTIN - 
ETAN * 

RHQRT3- 
END NAMELIST 
ENTER CHANGES 
•SINPUTS 

NAMELIST INPUT 

NSPDS - 15, 

APCNC ( I )» 


3.4000. 

4.2000, 


0 . 

0 . 

1 . 


1440, 


JCOOL - 0, 

0.049330, ETATTD- 

0.339500, XNRTD - 

0.487100, XMZXD - 

0. , TTIN - 

37.930000, FARGD - 

0.940000, R3QR2 - 

0 . 200000 , 

INPUT 

TO NAMELIST INPUT 


10 . 0000 , 
50.0000, 
*90 . 0000 , 
130.0000, 
* 20 , 

( I )« 

1 . 1000, 
1 .7000. 
2.4000, 
3.2000, 
4.0000, 
- 1 , 

( I >• 

1 . 0000 , 
0 . 

0.3300, 
* 1 . 


i 
5 
9 
13 
NPR 
APR 
1 
5 
9 
13 
17 
NAR 
ARN 
1 
5 
1 

NSTG 
DHQTD ■ 

TFFD - 
PSID - 
ANGSWX- 
PT1N ■ 

ETAN - 
RHQRT3* 

END NAMELIST 
NASA OUTPUT ON 


20 . 0000 , 

60.0000, 

100.0000, 

140.0000, 


. 2000 , 

.8000, 

.6000, 

4000, 


JCOOL * 
0.049330, 
0.339500, 
0.487100, 
0 . 

37.930000, 
0.940000 , 
0 . 200000 , 
INPUT 
IFC-15, 16 


4.2000, 


0. 

0 . 

1.1440, 

0, 


ETATTD 

XNRTD 

XMZXD 

TTIN 

FARGD 

R3QR2 


3 . 6000 . 3 . 

4.4000, 

0 , • • 

1.0000, 1- 


0.868000, 

1718 . 000000 , 

0 . 220000 , 

850.000000, 

0. 

0.291000. 


30.4404. 

40 

70.0400. 

80 

110.0000. 

150.0000. 

120 

1 . 4000 . 

1 

2 . 0000 . 

2 

2 . 8000 . 

3 

3 . 6000 . 

3 

4.4000. 

4 

0. 

0 

1 .0000. 

1 


0.868000, 

1718 . 000000 . 
0 . 220000 , 
850.0000001 
0 . 

0.291000, 


8800 , 
6000, 


0000 , 


. 0000 , 
. 0000 , 
. 0000 , 


.6000, 
. 2000 , 
.9000, 
. 8000 . 
.6000, 


. 0000 , 
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TURBINE FLOW 

FUNCTION 
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RPM, AND 
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7.0 ANALYTICAL BACKGROUKD 


The following section has been written in order to give the user a gen- 
eral idea of the type of turbine representation used in the program and the 
approach used in the derivation of the equations. Details of the deriva- 
tions together with sample calculations may be found in the Monthly Progress 
Reports (e.g., References 10, 11, and 12). 


7.1 TURBINE MAP REPRESENTATION 


Typically, cycle deck entry to a turbine map is through corrected speed, 
N/SQRT(T), and actual energy, DH/T, with turbine flow function, W*SQRT(T ) /P , 
and total-to-total efficiency being output. Total-to-total pressure ratio is 
sometimes used instead of actual energy as the second map entry. 

The discussion of the turbine map representation can be conveniently 
subdivided into two parts: the flow and the efficiency. 

The flow model is illustrated by the three sketches shown in Figure 4. 
The two curves on the top in the figure are used to generate the flow repre- 
sentation on the bottom. Sketch 4-1 shows the turbine stage characteristic 
(i.e., a plot of turbine loading against flow coefficient). Sketch 4-2 shows 
the dependence of the maximum value of the turbine flow function on corrected 
speed. With the corrected speed and DR/T known, the stage loading can be 
calculated, and the flow coefficient obtained from Sketch 4-1. Once the flow 
is choked (i.e., the choked branch of the stage characteristic), the flow 
coefficient remains constant for that speed. Sketch 4-2 is next used to ob- 
tain the maximum value of the turbine flow function at the corrected speed of 
interest. The value of the turbine flow function is then calculated. The 
equations used are as follows: 

U/A t - (2*R/60) (N//rRgg 0 T 

♦ - (DH/T)/(2(D//f) 2 /g 0 J] 

C/A* • (4»/Cos a 2 >/(U/A t ) 


1 
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Flow Coeff, $ ■ C 2 /U 
4-1 Stage Characteristic 


c 



4-2 Maximra Flow Function 



4-3 Flow Function Contours 


Figure 4. Turbine Flow Representation. 
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where 

DE/T • turbine stage specific enthalpy drop divided by ir.let total 

temperature, Btu/lbm R 

b/SQRT(T) ■ corrected speed, R?M/SQRT(R) 

R ■ pitch line radius, ft 

C ■ velocity, ft/sec 

• speed of sound at inlet total temperature, ft/sec 

TFF * inlet turbine flow function, lbm/sec*SQRT( R)/psia 

CosCoj) • nozzle exit angle 

The value of Cost^) i* obtained from the design point information as 
is the value of the pitch line radius. 

The efficiency model is illustrated by the five sketches shown in Figure 
5. Tne first three curves in the figure are used to generate the last two 
sketches. The "backbone" of the turbine map shown in Sketch 5-4 is the locus 
of the peak efficiency at each value of DH/T. This locus is obtained from 
Sketch 5-1, which shows turbine pitch line loading along the map "backbone" 
as a function of DE/T. The "backbone” efficiencies are obtained from Sketch 
5-2. This sketch gives the "backbone" loss (defined as the difference be- 
tween the ideal and actual values of DB/T) as a funeiton of DE/7. Uith the 
"backbone" loading and efficiencies known at each DK/T, Sketch 5-3 is used to 
evaluate the "of f-backbone" loss and to obtain the efficiency at any value 
of corrected speed. When the turbine "off-backbone" loss is plotted with the 
coordinates shown in Sketch 5-3, the resulting curves are nearly linear at 
any given DE/T. These five curves, three univariate and two bivariate, are 
sufficient to define the turbine map. Bote that as shown on Sketch 5-5, the 
design point does not generally fall on the "backbone" but is seperated from 
it by a "stand-off" distance which is calculated from design point informa- 
tion. 

The analytical basis for the five correlating curves is discussed in the 
following sections. 

7.2 TUK1IKE FLOW MODEL 

The turbine stage characteristic serves as the basis for modeling the 
turbine flow. An analytical expression for the stage characteristic of a 
turbine can be obtained by using the continuity, energy, and angular 
momentum equations, together with a number of relationships from the pitch 
line vector diagram. In deriving this equation, it is assumed that the pitch 
line flow angles at nozzle and bucket exit are invariant, and that the axial 
velocity ratio across the bucket is constant. 
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This later assuE 7 tlor, 1 e reasonable for incompressible or lov Mach number (M 0.25) 
flows. It should bt emphasized, however, that the final correlations have no such 
restrictions. 


The stage characteristic for a single-stage turbine with the above essurptio 
can be written in the form 


y = | ^ * 'L (x\) 


( 1 ) 


vhere 

4 i - DH/(2U2/g Q J) , Turbine stage Loading 

C, ■ Meridional velocity component 

* « C /U«, Flow coefficient 

22 2 

02 * Nozzle exit flow angle 

6 3 - Bucket exit flow angle 

R « Turbine Pitch Line Radius 


Selecting a reference point and assuming that 

c _5 *3 

tan a, 7 == T~ tan g 3 • a constant, 

*2 R 2 


this equatitin may be written as 



( 2 ) 


Tor mchoked flow, the refetenct point (indicated by the subscript 1 ) was 
selected st the design point on the turbine map. For choked flow, the refer- 
ence point was eelected st the critical point (subscript CS1T) on the ttfrbine 
map. The critical point is located at the value of Dfi/T at the reference 
apeed (e.g., 10C2) ^here nettle choking first occurs. For s variable-geometry 
turbine, the angle sj is calculated free the input value of nettle area ratio. 

The equation osed foT the normalised flow coefficient ia 

t C M fk t CF/^)g C/A t (W/^g m 

* (C 4 /A t ) t TJA (C/At)^ 1 1 A 

The velocity ratio, C/A^ for a point at a selected speed ia determined 
by calculating a pseud caret at rtiicb the Mach amber is aesmed to be mity 
(i.s., at the maximum turbine flow function for that speed). This pseudeerea 
it assumed te be constant for that speed. The Mach ember (velocity ratio) 
at any point on the speed line is then calculated from the usual flow function 
equations, by definition, the paeudnarea varies with speed in direct preper- 
tion te the maxism turbine flow fund ion variat ion with apeed. This method 
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of calculating the flow coefficient was found to give better results than 
that obtained using the first stage nozzle throat area for all speeds* For 
the case in which the first stage nozzle is choked, the two procedures are 
ident ical . 

Typically, the calculated and measured values of velocity ratio are 
within about 62, with the larger errors occurring at the higher hach numbers 
The predicted value of the normalized flow coefficient intercept (at DE/T*0) 
are within about 52 of experimentally derived values. These errors combine 
to yield flow errors on the order of 52 for nominal area maps (i.e., nozzle 
area ratios equal to 1). Slightly higher values may occur for other stator 
sett ings . 

A typical comparison between measured and calculated values of the tur- 
bine flow function is shown in Figure 6. This figure is for test turbine 
number 25 in Table 1. The maximum error shown on this figure is about 2.42 
and occurs at the lower end of the test data. This type of plot was used in 
obtaining the error estimates given above. 


7.3 TURBINE LOSS MODEL 


There are four key steps in the development of the equations governing 
the turbine off-design loss model. These steps are 

1. The development of an equation giving the turbine tot al-to-t otal 
efficiency at a general point in terms of nozzle and bucket effi- 
ciencies coupled with a semiempir ical loss term due to the depar- 
ture of the rotor incidence angle from the optimum. 

2. The transformation of the semiempir ical incidence angle loss law 
so as to eliminate the explicit occurrence of the incidence angle 
by introducing the stage loading. 

3. The differentiation of the resulting efficiency expression in 
order to obtain the locus of peak efficiencies. This peak effi- 
ciency ridge then becomes the "backbone 91 of the map. 

4. The substitution of the peak efficiency relationships back into 
the general efficiency equation in order to obtain an expression 
for the "off-backbone" loss. 

The development of the efficiency expression proceeds from the (h,s) 
diagram for adiabatic flow through a two-dimensional turbine stage as shown 
in the following sketch. 
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Using the station numbers shown in the above sketch, the definition of 
turbine total-to-total efficiency may be written in the form: 


t»tt 


AH 


AH + 


*o3 

*3 


(h 3 - h 3t ) 


♦ 



(h 3 - h2g) 


(4) 


By definition, the nozzle efficiency is equal to the ratio of the actuel 
nozzle exit kinetic energy to the ideal nozzle exit kinetic energy. A similar 
definition in terms of relative velocities holds for the bucket efficiency. 

For off-design calculations at any incidence angle (i), an additional loss 
term must be included, i.e., 

„ 3 . hj , . (j . J_) ^ (5) 

where 

* nozzle efficiency 
“ bucket efficiency 
W2 ■ bucket inlet reletive velocity 
W3 • bucket exit reletive velocity 
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The semiempirical incidence angle loss law is based on the assumption that 
the kinetic energy of the component of velocity normal to the optimum incidence 
angle is lost. This is a fairly standard assumption (see, for example, Reference 4). 
Powers other than 2 are frequently used on the cosine. 

Before introducing the incidence angle loss term into the efficiency expression, 
it was transfored into the following expression 

(/- O.s'U))^ - (W (. " 1 

'*<P i 5p- f 

The stage characteristic was then used to substitute for flow coefficient in 
terms of stage loading. The substitution of this results into the expression or 
efficiency yielded, after simplifying, the following results. 


( 6 ) 




J V\ 4 . vC 


K - 


i. L \ &''){%) esh 13 


(7) 


( 8 ) 


ft 


t„3 ¥> 


/z. 




U V 




(9) 


The temperature 'J “.^bET “ 

Sf to remain constant then Equation 7 can be differentiated, 

and thepeak efficienc y point loc ated^ 


4 


fu. 










■ i = 


L V<_ 


( 10 ) 


(ID 


then A=0 . 0 . 
expected. 


Note that if the noaale and bucket efficiencies equal unity in Equation (9) , 
1=0.0. Then there is no loss other than incidence and . ^ 0 p 
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By substituting Equation* 10 and 11 into Equation 7, the following ex- 
pression for "off-backbone 11 loss can be obtained. 

<i2 

Equations 10 f 11, and 12 give the location of the peak efficiency, the 
magnitude of the peak efficiency, and the variation in efficiency as vc move 
away from the peak. 

These equations represent the stage loss characteristic of a turbine. 
The design point information is used to obtain the initial values of A and B 
as well as the values of the blade row efficiencies and metal angles. The 
loss equations are then applied at incremental values of DHQT starting at 
aero to obtain the turbine efficiencies. Approximate relationships are used 
for temperature and velocity ratios to obtain new values of A and B for each 
DHQT. 


Typically, the calculated values of the loss slopes and those obtained 
from air turbine test data are within about 5Z for corrected speeds within 
plus or minus 20Z of the design point value. The values of the "backbone 
efficiencies" generated by the above equations do not include either Reynold's 
Number effects or the severe rotor exit losses encountered near exit annulus 
choke. In order to account for these effects, the loss along the peak effi- 
ciency ridge was empirically modified using the results of the NASA air 
turbine tests. Although relatively good correlation existed between the 
different test turbines at low values of DHQT, the drop in efficiency in the 
neighborhood of exit annulus choke was so severe thst correlstion was diffi- 
cult. For this reason, variations in efficiency on the order of SZ can be 
obtained in this region. 

A comparison between measured and calculated values of the turbine total- 
to-total efficiency is shown in Figure 7. This figure is for test turbine 
number 25 in Table I. The maximum error shown on this figure is about 0.8Z. 
This type of plot was used in obtaining the error estimates given above. 
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7.4 Comparison of Radial Turbine Test Results and Program Output 


A comparison was made between the program output and the "SOLAR' 1 radial 
turbine described in Reference 20. The key design point parameters are summarized 
in Table IV. These values are used as the default settings for radial turbines as 
shown in the second example in Section 6. 

The performance map is shown in Figure 8. The exit annulus choke curve forms 
the physical limit of the map. The portion of the map above this limit line re- 
presents a mathematical extension of the map which is necessary for cycle deck 
iteration purposes. The predicted swirl map is shown in Figure 9. The predicted 
equivalent weight flows are compared with the test data in Figure 10. The solid lines 
are the predicted values. In Figure 11 the total-to-static efficiency is compared 
with the so-called blade to jet speed ratio. This parameter is the ratio of the in- 
let wheel speed to the velocity calculated by expanding the flow isentropically 
through the inlet total to exit static pressure ratio. This is the usual manner of 
presenting radial-turbine efficiency data. The solid lines are the predicted values. 
The sharp breaks in the lines are due to interpolation. Values of DH/TA at intervals 
of 0.005 were used to read the map efficiencies and flows from Figure 8. The breaks 
are due to missing the peak efficiency at a given speed. A smaller interval could not 
be used due to size limitations on the program. 

In general the comparison is quite good, and is probably within the accuracy 
that the data could be read from the relatively small curve given in Reference 20. 
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TABLE VI . 

Design Point Data - INFLOW Radial Turbine 
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SUIRL MAP 
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blaoe-jet speed ratio 

(Solar 4.75-in. Dia. Radial Turbine) 


Figure 11. Radial Turbine Efficiency Comparison. 
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LIST OF SYMBOLS 


A, a 

B, b 
C 


Cp 

to 

h 

i 

J 

N 

P 

R 


S 

T 

U 

V 


Const ant • 

Constants 
Velocity, fps 

Constant Pressure Specific Heat, Btu/lbo*R 
Diaens iooal Constant, 32,17 ft lbm/lbf sec 2 
Enthalpy, Btu/lbo 

Incidence Angle (i * * ®2op^ • degrees 

Mechanical Equivalent of Heat, 778, 16 ft. lbf/Btu 

Speed, rpo 
Pressure^ psia 

Radius, ft 

Gas Constant, 53. 35 ft lbf/lba # R 

Entropy, Btu/lbm"R 
Absolute Temperature, # R 
Wheel Speed, fps 
Relative Velocity, fps 


a 

6 

AH 

r 

x 

0 

6 

* 

n 


Angle of Absolute Velocity With Axial, degrees 

Angle of Relative Velocity Vector With Axial, degrees 

Prop in Total Enthalpy, Btu/lbn 

Ratio of Specific Beats 

Enthalpy Loss Coefficient 

Fluid Density, Ibm/ft^ 

Ratio or Total Temperature to Standard Temperature 
Turbine Stage Loading {P ■ AH/(2ti 2 /g 0 J)] 

Flow Coefficient (♦ " C r /U) 

Ef f ic iency 


Subscripts 

1 

2 

3 

B 

d 

N 

0 

op 

P* 

rr 

s 

s 


Nozzle Inlet 
Rotor Inlet 
Rotor Exit 
Bucket 

Design Point Value 
Nozzle 
St agnat ion 
Opt iaua 

Peak 

Total-to-Total 

Isentropic 

Axial 
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